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DISCOVERING  HIDDEN  TOTALLY  LEONTIEF  SUBSTITUTION  SYSTEMS 

by 

George  B.  Dantzlg  and  Arthur  F.  Velnott,  Jr. 

Abstract 

A constructive  procedure  is  given  for  determining  the  existence 
of  and  evaluating  (when  it  does  exist)  a nonsingular  matrix  that  trans- 
forms a system  of  linear  equations  In  nonnegative  variables  Into  a 
totally  Leontlef  substitution  system.  The  computational  effort  in- 
volved Is  about  that  required  to  optimize  the  given  m-row  linear  system 
with  m+1  different  linear  objective  functions. 


DISCOVERING  HIDDEN  TOTALLY  LEONTIEF  SUBSTITUTION  SYSTEMS 


George  B.  Dantzig  and  Arthur  F.  Veinott,  Jr. 


The  system  of  m linear  equations  in  n nonnegative  variables 


Ax  = b , X > 0 


Is  called  a Leontief  substitution  system  [2]  if  (i)  each  column  Aj 
of  A has  at  most  one  positive  element,  (11)  b » 0 and  (111)  the 
set  of  solutions  to  (1)  Is  nonempty.  If  also  that  set  Is  bounded, 

(1)  Is  called  a totally  Leontief  substitution  system  [6] . In  either 
case.  It  Is  known  that  A has  rank  m.  Such  systems  are  discussed 
in  [1]  - [8]. 

Saigal  [8],  [7]  calls  (1)  a hidden  totally  Leontief  substitution 
system  if  there  exists  a nonsingular  matrix  n such  that 


(IIA)x  = rib  , X > 0 


is  a totally  Leontief  substitution  system.  The  purpose  of  this  paper 
Is  to  give  a constructive  method  for  determining  whether  or  not  (1) 
has  this  property,  and  if  so,  to  find  11  . 


Substitution  Classes.  Associated  with  any  feasible  m x m basis 
® ^ “ l,...,m  , the  set  of  column  indices 

J such  that  Aj  , if  substituted  for  B^  , forms  a feasible  basis. 

Each  substitution  class  is  nonempty  since  it  Includes  a j with 


1 


Aj  ■ . In  general  the  depend  on  B and  b and  can  be  overlap- 

ping when  there  are  degenerate  basic  feasible  solutions.  For  the  totally 
Leontief  substitution  case,  however,  the  are  independent  of  the 

choice  of  B and  b >>  0 ; indeed,  consists  of  all  j such  that 

Aj  has  a positive  element  in  the  same  row  as  B^  • Also  the  par- 

tition the  column  indices  l,...,n  . Thus  every  submatrix  B consist- 
ing of  m columns  of  A with  a positive  element  in  each  row  forms  a 
(nondegenerate)  feasible  basis,  and  conversely. 


The  Algorithm. 

Step  1.  Find  a feasible  m x m basis  B and  determine  substitution 

classes  S, ,...,S  with  respect  to  B , b . Terminate  if  there  is  no 
1 m 

feasible  basis  or  if  the  substitution  classes  do  not  partition  the  col- 
umn indices  of  A . Otherwise  go  to  Step  2. 

Step  2.  Solve  the  linear  program  of  maximizing  z = Z x^  subject  to  (1), 
Terminate  if  z is  unbounded  above.  Otherwise  go  to  Step  3. 

Step  3.  For  each  i = l,...,m  , determine  the  i^^  row  of  11  as  any 
vector  such  chat 


(3) 


TT^b  > 0 


tt^A^  1 0 for  all  j ^ 


Terminate  if  for  any  i •=  l,...,m  the  system  (3)  is  infeasible.  Other- 
wise terminate  with  11  . 
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Theorem.  If  the  algorithm  terminates  with  11  , then  IT  is  nonsingular 
and  (2)  is  a totally  Leontief  substitution  system.  Otherwise  (1)  is 
not  a hidden  totally  Leontief  substitution  system. 

Proof.  If  (1)  is  a hidden  totally  Leontief  substitution  system,  then 
there  is  a feasible  basis,  the  associated  substitution  classes  partition 
the  column  indices  of  A and  z in  Step  2 is  bounded  above,  because 
these  properties  are  invariant  under  nonsingular  transformations  IT  . 
Also  there  is  a nonsingular  matrix  IT  whose  i^^  row  satisfies  (3) 

for  each  i . Thus  if  the  algorithm  terminates  without  obtaining  II  , 
then  (1)  is  not  a hidden  totally  Leontief  substitution  system. 

If  the  algorithm  does  terminate  with  IT,  then  HA  has  at  most 
one  positive  element  in  each  column  (from  Steps  1 and  3) , lib  >>  0 
(from  Step  3)  and  (2)  has  a solution  (from  Step  1),  so  (2)  is  a Leontief 
substitution  system.  Hence  IIA  has  rank  m,  implying  IT  is  non- 
singular  and  so  (1)  and  (2)  have  the  same  solution  set.  Thus  the 
boundedness  of  the  solution  set  of  (1)  implies  that  is  so  of  (2) , so 
(2)  is  a tot  illy  Leontief  substitution  system. 


Computational  Remarks.  The  computational  effort  required  to  execute  the 
algorithm  is  about  that  required  to  solve  the  linear  program  of  minimiz- 
ing cx  subject  to  (1)  with  m+1  different  objective-function  vectors 
c “ determine  the  substitution  classes  in  Step  1 requires 

computing  b'  • = B ^b  and  A^  • ^^jk^  " ® ^j  each  j . Then 

the  substitution  classes  partition  the  column  indices  of  A if  and 
only  if  for  each  j there  is  a unique  k ■ k(j)  that  minimizes 
subject  to  Ajj^  > 0 . In  that  event  J ^ ^lt(j)  each  j . Step  2 
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involves  solving  the  linear  program  with  « -1  for  all  j . Finally, 
Step  3 necessitates  solving  m linear  programs.  The  i^^,  1 < 1 < m , 
of  these  has  * 1 for  all  3 ^ ^ otherwise.  If  optimal 

simplex  multipliers  exist  therefor,  they  satisfy  (3).  If  no  such 

multipliers  exist,  (3)  is  infeasible.  Incidentally,  Step  3 can  be  stream- 
lined somewhat  by  modifying  the  i^^  linear  program  so  that  all  but  an 


(arbitrary)  one  of  the  variables  x^  with  3 ^ omitted. 


A 
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